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We present results for the spectrum of bb bound states in the quenched approximation for three different values 
of the lattice spacing. Results for spin-independent splittings are shown to have good scaling behaviour; spin- 
dependent splittings are more sensitive to discretisation effects. We discuss what needs to be done to match the 
experimental spectrum. 



1. INTRODUCTION 

Accurate calculations of the hadron spectrum 
in lattice QCD require control of systematic er- 
rors, and a major source of these errors arises 
from the use of a space-time lattice with finite lat- 
tice spacing. Recent understanding of these dis- 
cretisation errors has meant that the 'brute force' 
approach of extrapolating to zero lattice spacing 
has been replaced by the requirement of improv- 
ing the lattice action to the point where results 
at different values of a agree to within specified 
precision. This also allows us to use effective ac- 
tions for QCD for which a continuum extrapola- 
tion cannot in any case be done, but which are 
nevertheless good for particular physics. 

The approach of Non-relativistic QCD fits into 
this picture ||-||. The NRQCD action is well- 
suited to a description of non-relativistic quarks 
such as b quarks within bottomonium or Wight 
mesons In fact we believe that it is the best 
approach for these systems. The operators of 
the action are classified by the powers of v 2 /c 2 
that they contain with coefficients that can be 
matched to full QCD in perturbation theory. The 
renormalisability of QCD is lost, however, so an 
explicit momentum cut-off (non-zero lattice spac- 
ing) is necessary for NRQCD and the coefficients 
will contain radiative corrections which diverge as 
powers of this cut-off, A/itiq, so A/ttiq — > oo is 
not possible. Within a range of A around run 
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physical results independent of A are possible, 
however. Here we will discuss to what extent this 
situation is realised for current results. 

1.1. Simulation Details 

We use an evolution equation for quark propa- 
gators 
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Ho is the lowest order (in v 2 /c 2 ) term in the 
Hamiltonian, the kinetic energy operator: 
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The correction terms that we include in SH are 
0(w 4 /c 4 ) ||. They comprise relativistic correc- 
tions to the spin-independent Hq as well as the 
first spin-dependent terms that give rise to spin- 
splittings in the spectrum. 
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The last two terms in SH come from finite lattice 
spacing corrections to the lattice laplacian and 
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the lattice time derivative respectively Q2J. We 
tadpole-improve our lattice action by dividing all 
the Us that appear by Uq, which we take from the 
fourth root of the plaquette (u p). We then work 
with tree-level values for the a, i.e. 1. 

Table 1 shows the parameters used in the calcu- 
lation at 3 different values of /3, for the standard 
Wilson gauge action ||. 
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Table 1 

The parameters used. We are grateful to the 
UKQCD collaboration and to Kogut et al for the 
use of their gauge field configurations. 



2. RESULTS 

2.1. Radial and orbital splittings 

Radial and orbital splittings (when spin- 
averaged to remove spin effects) are calculated 
at next-to-leading order both in terms of a non- 
relativistic expansion and in terms of discreti- 
sation corrections by the action of eq. [| We 
find very little remaining a dependence when 
we take dimensionless ratios of spin-independent 
splittings, as in Figure |l|. The disagreement with 
the experimental results, shown as lines is pre- 
sumably an error from the quenched approxima- 
tion since higher order (physical) relativistic cor- 
rections should have a ~ 1% effect. This inde- 
pendence from the lattice spacing allows us then 
to extrapolate quenched and dynamical (rif — 2) 
results to the real world 

2.2. Spin splittings 

The spin splitting that can be calculated most 
accurately is the hyperfine splitting between T 
and ?7f,. No experimental value is known and so 
the lattice results can make a useful prediction if 
we can reduce systematic errors below 10%. Our 
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Figure 1. Various spin- independent ratios of 
splittings as a function of lattice spacing, set by 
the Xb — T splitting for NRQCD in the quenched 
approximation. For A(25— 15)/A(1P— IS) two 
values are shown. The squares are the raw results, 
the diamonds include 0(a 2 ) corections for errors 
in the Wilson plaquette gauge action, estimated 
from the hyperfine splitting 0. 



aim here is to check the scaling of the quenched 
results before attempting an extrapolation in n/ 

0- 

Because the spin splittings are sensitive to the 
heavy quark mass it is important that this is 
tuned correctly. For our simulation parameters 
the quark masses we use are closest to the cor- 
rect ones if we fix the scale from the 2S-1S split- 
ting. Figure || then shows the hyperfine split- 
ting in MeV using this scale at the 3 different lat- 
tice spacings. Our results are given by the plain 
squares, and indicate strong scaling violations, as 
might be expected of such a short distance quan- 
tity with 0(a 2 ) errors. The fancy squares are 
rescaled by (uqp/uol) 6 to give an indication of 
the results that would be obtained by using uq 
from the Landau gauge link rather than the pla- 
quette. They are rather flatter. 

The plain diamonds show results from Manke 
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Figure 2. The hyperfine splitting in physical units 
vs a 2 (setting the scale from the 2S-1S splitting) 
for various NRQCD actions in the quenched ap- 
proximation. See text for explanation of the sym- 
bols. 



et al j8[ using an action improved to the next 
order for spin-dependent terms This means 
adding additional relativistic and discretisation 
corrections which however will not change the 
spin-independent spectrum (or lattice spacing). 
Since these calculations use the same parameters 
as ours, we have shown them in Figure using 
our 2S-1S splitting to set the scale. The plain di- 
amonds show much better scaling behaviour than 
the plain squares and illustrate the fact that the 
discretisation and relativistic corrections act in 
different directions The discretisation cor- 
rections tend to increase the hyperfine splitting 
(obvious from the behaviour of the plain squares) 
but the relativistic corrections tend to reduce the 
splitting (so that the plain diamonds are flatter 
and below the plain squares). 

The fancy diamonds illustrate approximately 
what would happen to the plain diamonds if uql 
had been used instead of uqp, using the leading 
order rescaling above. Now very good scaling is 
observed and this is encouraging, but must be 



checked in a complete calculation using uql ■ Note 
that the SESAM collaboration have used exactly 
this action with extra spin-dependent corrections 
and uql in their comparison of results with n/ = 
and 2 JK]]. 

Statistical errors in the p fine structure in cur- 
rent calculations make it hard to reach clear con- 
clusions about the size of scaling violations ||. 

3. CONCLUSIONS 

NRQCD actions can give physical results for 
bottomonium splittings which scale accurately for 
a reasonable range of lattice spacings. This is 
true for spin-independent splittings for the action 
of equation ^; for spin-dependent splittings it is 
likely that this will happen with the next order 
of relativistic and discretisation corrections |@,ffol . 
At this level radiative corrections to the leading 
order terms also have an effect and it will be nec- 
essary to decide which uo to use and what the ra- 
diative corrections beyond tadpole-improvement 
are Q. 
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